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1 Introduction 

The de Sitter spacetime is the (one and only) positively curved maximal 
symmetric spacetime, in other words the positively curved spacetime that has 
the same number of Killing vectors as Minkowski spacetime. As such, is an 
ideal candidate for developing a quantum field theory on a curved but fixed 
background, covered with a certain coordinate chart. Following the steps 
of the canonical quantisation of the flat case, there has been constructed a 
quantum field theory on the spatially flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) chart. [1] 

The Killing vectors of the manifold give rise to the generators of the 
representations of the isometry group carried by spaces of matter fields. |2] 
Therefore, the generators are operators related to the conserved quantities 
given by the Noether theorem, and from them one constructs the conserved 
one-particle operators of the corresponding quantum field theory. 

The generators also commute with the operator of the field equation. In 
this respect, they or any other operator which commutes with that of the 
field equation are called conserved operators. It is important to point out 
that a set of mode functions of a specified field equation of given spin on 
a spacetime (covered with a certain chart) is completely determined by the 
complete set of commutative operators (CSCO) {S, A, B, C} chosen to label 
the modes. These must include the operator of the field equation S, and 
some operators A,B,C, of which the modes are eigenf unctions, the number 
of which depends on the spacetime dimension and if the fields have spin. For 
scalar fields on a 3+1 dimensional spacetime, 3 such operators are needed. 
If these operator are differential, then their eigenvalues can be thought of as 
separation constants for the Klein-Gordon equation. For fields with spin, an 
additional operator is required, which encodes in some way the spin. |3] 

Though the scalar field is widely used in cosmology |H [5] , there hasn't 
yet been found evidence of a scalar elementary particle. Nevertheless, the 
free scalar field can be a good model for elucidating the inner workings of a 
field theory. E] 

A number of solutions for the Klein-Gordon equation in different charts 
on the de Sitter spacetime have been found, such as the static chart |9j, the 
hyperspherical one [10], or the spatially- fiat FLRW one. [11] These three 
solutions are the best known, being treated even in classical texbooks. [12] 

There have been computed even Bogolyubov coefficients between the 
modes on the above mentioned charts. ^13j Moreover, particle creation can 



be explained as emerging from non-null /3-Bogolyubov coefficients, especially 
due to the thermal nature of the de Sitter horizon, and a number of studies 
have been conducted in this regard. |T31 115] 

Other new modes can be defined either in new charts, either with new 
operators on the old charts. Recently, the theory of free scalar fields on de 
Sitter spacetime has also been discussed in Beltrami charts. [K] And a new 
hamiltonian-like operator has been on the spatially fiat FLRW charts. pQ 
In fact, two spatially fiat FLRW coordinate charts are needed in order to 
cover the whole of de Sitter space. These two charts don't overlap, one of 
them describes an expanding spacetime to which we refer as "the de Sit- 
ter expanding universe", and the other a contracting one. The mode solu- 
tions on the spatially-flat FLRW chart are also momentum (or momentum 
squared) eingenvectors (called "momentum basis modes", either plane- wave 
modes- eigenfunctions of {SKGyP}Ul\y *^^ spherical ones- eigenf unctions of 
{Skg, P^, L"^, Lz}, as used in [13]). 

Recently, Cotaescu et. al pQ managed to find, through a technique involv- 
ing time evolution pictures on spatially fiat FLRW manifolds, scalar plane 
wave solutions that are eigenfunctions of a newly-defined conserved operator, 
a modified Hamiltonian H (modes of which are called "energy basis" plane 
wave modes, and are eigenfunctions of {Skg, H, Mi}). 

They also argued that, due to the uncertainty between the energy and 
momentum on de Sitter spacetime, both momentum basis and energy basis 
quantum modes are required in order to build a complete de Sitter quantum 
field theory, for which some applications have already been worked out.|17[ 
H] 

Still, the spherical energy basis scalar quantum modes (eigenfunctions 
of {Skg, H, L^, L^}) remained unsolved. This paper is devoted to finding 
the analytical form of these scalar modes, on the spatially flat de Sitter 
FLRW chart, via the Schrodinger picture (SP) of the time evolution technique 
developed specially for spatially flat FLRW spacetimes. [19j 

In the second section, we make a brief review of the known results, with 
the formalism used afterwards to derive the new results. This is done not just 
to accommodate the reader with the already known solutions or used conven- 
tions, but is necessary in order to understand both the similarities and mainly 
the differences between the momentum basis quantum modes and the cor- 
responding energy basis ones. The second section concludes with expressing 
the planar energy basis modes in a new form, as a sum of two hypergeometric 



functions. The third section contains entirely original work. It consists of the 
first derivation of the new spherical scalar energy-basis modes, and a way of 
determining them through a different method of the one used to obtain the 
plane-wave correspondents. [Ij The paper concludes with the fourth section 
where the transition coefficients between all the aforementioned modes are 
briefly computed, showing that the energy basis modes are indeed different 
from the previously known momentum basis ones. 

2 Previous solutions on the spatially fiat FLRW 
chart 

2.1 Momentum basis plane waves 

The Klein-Gordon equation, on the 3+1 dimensional de Sitter spacetime, 
in the spatially flat FLRW chart with line element ds^ = dt^ — e^'^^dx ■ dx, 
parametrised by (spatial) Cartesian coordinates {t, x} takes the form 

{d^ + 3ujdt - e-'"*A,,j,,, + m^) $(t, f) = 0. (1) 

The first quantum modes ever written on the de Sitter spacetime, were 
the eigenf unctions of the momentum operator, in this chart. [TT] By promptly 
solving the three eigenvalue equations Pi^(t, x) = Pi^(t, x), one obtains that 
the fundamental solution of the Klein-Gordon equation is in this case sepa- 
rable as f^(t,x) = N fp(t)e^^'^ , which leads to the equation for the temporal 
part: 

^/,(t) + 3u;|/p(t) + ip'e-'-' + m')U{t) = 0, (2) 



which has the solution: 



Uit) = Ce'-^Hi'^ (^e--) . (3) 



with z/ = i^/i2 - 9/4. 

It is worth pointing out that the choice of the function Hu as the solution 
of the Bessel equation is due to the fact that one chooses to work with 

(2) 

progressive waves, rather than regressive ones {Hi - corresponding to /*). 

The fundamental solutions of the equation can be normalisation with the 
help of the scalar product expressed in the Cartesian chart as 



{f{t,x)J'{t,x)) = ie'-' d'xnt,x)dtf'{t,x), (4) 

Jrs. 

which involves the Wronskian of two Hankel functions 



'^ \uj yy ' \ui y n 



Here, two cases are included: when u is real and fi ^ 3/2, and when u is 
pure imaginary and fi > 3/2, where an extra factor e~^'"'^ appears, being 
just the Heaviside step function used to encompass both cases. With these 
in mind, the normalised modes have the following expression: 

2.2 Momentum basis spherical waves 

The Klein-Gordon equation on this chart ([1]) can be put in spherical coor- 
dinates. The momentum basis for spherical waves is given by the CSCO 
{Skg, P^i L'^1 Lz}, so its fundamental solution must be of the form 

/p,,^,(t,r,^,0) = iV/p(t)i?,,,Wn-,(^,0), (7) 

where the radial part is given by the eigenvalue equation P'^Rp^iir) = p^Rp^iir). 
Using the Laplacian in spherical coordinates leads to a spherical Bessel dif- 
ferential equation: 

^2RvAr) + -^^R,Ar) - [^-r^ - P') RvAr) = 0. (8) 

The solution can be written as Rp^r) = Cji{rp) = C"4^Jj_,_i (rp), where the 
temporal part fp{t) is the same as for the plane waves ([3]), such that the 
normalised modes are in this case 

A y CO y/r \uj y 

(9) 
by using the scalar product expressed in spherical coordinates, which reads: 



_r>00 



{f{t,r,e,<f)),f'{t,r,e,<P))=te'''' / drr' / dQr{t,r,e,<P)dtf{t,r,e,<f)). 

Jo JS2 

(10) 



It should be noted that the spherical moinentuin modes can also be ob- 
tained from the plane wave ones, by performing a Rayleigh expansion- ex- 
panding the plane wave into a sum of spherical waves: 



^^q■x 



^ CO / 



qr 



2.3 Energy basis plane waves 

Cotaescu et al. found the scalar modes that are eigenfunctions of a new 
Hamiltonian. [1] Since the —idt is not conserved in the sense that is doesn't 
commute with the field equation operator, the energy basis is defined as 
being given by the CSCO {Skg: H,Mi}, where H = —idt + ujx^di is a modi- 
fied conserved Hamiltonian, and A/i are any two of the three non-differential 
operators which encode the momentum vector's direction. 

Passing from the natural picture (NP) to the Schrodinger picture (SP), 
according to the rigorous procedure fully described in pQ, the Klein-Gordon 
equation of the natural picture ([1]) becomes in the Schrodinger picture, in 
the same chart: 

[{dt + cox'dif + 3u{dt + ux%) - A,,j,,, + m^] ^^(t, f) = 0. (12) 
Expanding the field as $5(x) = $^ (x) + $^ (x) with 

/CO P 

dE / £q ^^^\e, g)e^*(^*-^^-^), (13) 

one arrives to the Klein-Gordon equation in momentum representation, in 
the SP: 

{±zE + u{q%^ + 3))' - 3u {±iE + u{q'dq^ + 3)) + g^ + m'] ^f\E, g) = 0, 

(14) 
where for example for the positive frequency part $5. (£", q) = hs{E, q) a{E, n), 
and the function hs being of radial type, it satisfies the equation 

[{±tE + ooiqdg + 3)f - Sou {±iE + uj{qdg + 3)) + q^ + m"] hf\E, q) = 0, 

as) 

which can be shown to have solutions of the form: /^^(e, s) = Cs~'^'^~^^'^Hi, (s), 
where the variables e = E/u,s = q/u and /i = m/u have been rescaled for 
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simplicity. Again, the choice Hi) as the solution of the Bessel equation is 
made due to the convention that we express / using progressive waves. 

From the desired expansion of the scalar field in SP, the expression for 
the mode functions (in SP) can be read off as: 

fl.{t,x)=Ne-^^' dsV^ss-^^Hi'\s), (16) 

Jo 

and by passing back to the one in NP, followed by normalization the energy 
basis modes (in integral representation) are obtained: 

(17) 
This is the form in which the plane-wave energy basis scalar modes were 
first given. [T] The authors searched for the plane waves which are eigenfunc- 
tions of a CSCO that included non-differential operators Mi. Due to this, one 
could not simply solve the equation by using the extra eigenvalue equations, 
and had to resort to writing the equation in the momentum representation. 
However, as we shall see, this is not the case for the energy basis spherical 
waves. 

Furthermore, the integral in the expression of the modes can be solved, 
giving a Gauss hypergeometric function. 



fE,n{t, x) = ^e 2 e 2 



Ljt\U 



w„-%t^isie(u-3/2) {—2ie "^ ) 



-Ult\ 



11 (^— iwn • a; — ze~""'-)~*^+^/^+'^ 

In this paper, we present a new explicit form for the above expression 
of the energy basis scalar modes, which can be written in terms of other 
hypergeometric functions [21] as: 

r(a+)r(a_)2Fi (a^, a_; i; co\n ■ xfe'^' 



UJ ^E^_iEt 



fEAt,^) = j^:;^e 2.e 



-2un ■ xe'^*r ( a+ + ^ V L + i") ^F, (a+ + i, a_ + i; | co\n ■ xfe^^' 



(19) 



where a± = | — ^ ± |. This is the form which can be compared with the 
new results that will be deduced in the next section. 



3 Energy basis spherical waves on the spa- 
tially flat FLRW chart 

We now proceed to the main original part of this work- finding the spherical 
energy basis quantum modes. These are defined as eigenfunctions of the 
CSCO {£kg^H,U',Lz}. Unlike the plane wave energy basis modes, the 
spherical ones being eigenfunctions of differential operators, can be found 
by solving the Klein-Gordon equation, in the spherical FLRW chart, using 
the associated eigenvalue equations. One of the main points of the paper is 
showing how these modes can be found, by solving the equation in this way. 
In the Klein-Gordon equation in SP (I12p . we pass from the Cartesian 
coordinates {t, x} to spherical coordinates {t, r, 6, 0}, so the equation changes 
to: 

{dt + urdrf + ^u{dt + urdr) - dl --dr-^+ nA ^s{t, r, 9, 0) = 0. 

(20) 
The "energy basis spherical" quantum modes we wish to find are solutions 
of this equation, but (as their name implies) they are also eigenfunctions 
of the following operators: H,L^,Lz. By promptly solving the eigenvalues 
equations, while being in the SP, for the CSCO H = —idt, L^ = — Ag <^, L^ = 
— A<^, it is straightforward to show that the solution of the equation must 
include the factor e~*'^*, and the spherical harmonics Yj,m,(6',0), such that 
the fundamental solutions are separable as: 

fE,i,m,it,r,e,<l)) = iVe-^^*i?|,,(r)l^,^,(^,0). (21) 

It is of note that the separation of variables occurs only here in the SP, and 
only thanks to the use of this time evolution technique. [T] We stress again 
that in SP, H = —idt is conserved, while in the NP only the modified Hamil- 
tionian H = —idt — oox^di is conserved. By inputting the above solution, and 
making use of the chosen operator eigenvalue equations, the Klein-Gordon 
equation becomes an ordinary differential equation for the radial part: 

(cjV^ - 1)— - + [Au^r - 2iuEr - - W 
dr^ \ r dr 



-E^-3zcoE+^-^^^ + m' 



R'{r) = 0. (22) 



Changing the variable from r to z = u'^r^, and dividing by ~4u'^, the 
equation becomes: 

2 Q- in I 1^ ,,2^ 



d'^ f 5z 3 \ d fe^ 3ze /(/ + 1) n' 



RU^) = 0, 



(23) 
where for brevity e = E/uj and /i = tti/w. 

Making an ansatz for Rs{z) as R^ i{z) = z^S^^^z), we obtain an equation 
which has the form of the Gauss hypergeometric equation 



(24) 






/2 3/ 3ze ze/ e^ /i2 

+ '4-4+T + T+4-T 



sU^) = 0. 



Discarding the solution that is singular at the origin, we remain only with: 



/ 



SeAz) = 2F1 cr+ - -, o-_ - -; / + -; 2; 



(25) 



where again (j± = | — ^ ± | and v = "^-v/yU^ — f • 

Bearing in mind all of the above, and by passing back from the variable 

z to r: 



(26) 



fl,^,it,r,e,<l>) = Ne-^'^\ujry,F,{a 



and then from SP to NP, we obtain the following solution to the Klein-Gordon 
equation: 

/ , 3 ,.^.^_, , 

(27) 
Next, these quantum modes will be normalised. In order to do that, part 
of the solution must be written as a Hankel transform of a certain function: 

1 



fE,i,m,{t,r,e,<P) = iVe-*^*(c^re-*)'2i^i (cj+ - ^,a_ 



-;/ + -;a;VV-lF„ 



fE,i,n.,{t,r,e,(j)) = Ne-'^\ioe''y^Yi,^X^,<p)Hx[9{oor)], (28) 



where 



UM^r)] = r^n.F, ( a+ - ^, a_ - ^; / + I c^V^e^^* 



J+h 



(29) 



According to [21] , this is the Hankel transform of order A = / + | of : 



g{s) = u^ 



which can be evaluated as 



X- 



j^s-''-^Hl^\se-'''), (31) 



r(a+-i)r(a_-|)' 

such that 

As a result of this, the quantum modes in integral representation can be 
written as: 

r (i -L I'] 



X 



rK-i)r(a_-i) 

F,,„,(^,0)^ / s-''Hi^\se-^')Ji^i{ujrs)ds. (33) 



Applying the scalar product, written in spherical coordinates (TTOj) . and 
taking note of the closure relations and the Wronskian of the Hankel functions 
(ISl), the normalisation constant is: 

ef r ((T+ - ^) r ((T_ - ^) , , 

N = —— — ^^-- j^ — ^v ^, 34 

such that the modes have the simple expression in the integral form (dropping 
the phase factors): 

fE,i,mAt,r,e,<l>) = -i=e-^e^®(^-t)ll,^,(^,0) 



X 



-^j^ s-Hl'\se-')J,^.{u;rs)ds, (35) 
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or the full, normalized explicit expression 

2\/27TUJ2 i [I' + 2) 

X2F1 (^+ - ^' ^- - ^; ^ + I ^'r'e'^') yi,mX()^ <P)- (36) 

The completeness relation for these mode functions can be verified with 
their integral form. 

Alternatively, one can use the Rayleigh expansion f lTT|) on the plane wave 
energy basis normalised modes f lT7|) to obtain the corresponding spherical 
ones ( 15^ . which ensures the agreement between the above approach and the 
method of [T]. 

These modes are important because they can be used as scalar solutions 
of given energy in scattering problems like partial wave analysis on de Sitter 
spacetime, where for instance an incoming plane wave of either definite mo- 
mentum dnj or definite energy of the form (IT7|) scatters off a central potential, 
resulting in a spherical wave of definite energy f l35p . This could shed light 
on the problem of correctly defining the scattering cross section on de Sitter 
spacetime. 



4 Transition coefficients 

One can compute the scalar products between any of the quantum modes 
fl6|) ■ fl9|) . flT7|) and fl35|) . For the scalar product between plane and spherical 
waves, it is necessary to expand the plane ones into spherical waves, via the 
Rayleigh expansion formula, in order to be able to use the scalar product in 
spherical coordinates (ITU]) . 
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These transition coefficients have the following expressions: 



3 



{fy, fE,n) = -^S'i^ - ^p)e-^" ^^- , (37) 

{fp, fp,l,m,) = -5ip - \p\)H)%rr,iO,^, <Pp) , (38) 

_3 
{UJE,l,m,) = -^{-t)%m,ie,^Ap)e''^'''^, (39) 

_ 1 
{fE,nJp,Wm,) = -^^{-l)%,raXOHAn)e'^'''-, (40) 

UE,njE',l,m,) = 5{E - E'){-l)%^Xdn.(t>n). (41) 

{fp,Lmi, fE,l',m') = r- ^l.V^miM^'^' "^ " , (42) 

V27ra; 

(43) 

while their converses satisfy: 

(/2,/i) = (/i,/2r. (44) 

All these transition coefficients can be thought of as a-Bogolyubov coeffi- 
cients. However, the /3-Bogolyubov coefficients between any of the discussed 
mode functions cancel out, thanks to the fact that the Wronskian of two 
Hankel functions of the same kind is null. 

(/i,/2) = 0. (45) 

This means there is no mixing of positive and negative frequencies. 

5 Conclusions 

Scalar spherical quantum modes on the de Sitter expanding universe, that 
are eigenf unctions of a modified Hamiltonian operator have been computed, 
completing the free scalar field mode analysis on de Sitter spacetime in the 
spatially flat FLRW chart. [1| Our scalar modes, in their hypergeo metric 
form remind of the ones found long ago on the static chart. [9j An explicit 
alternate form of the plane-wave energy basis quantum modes, in terms of 
two hypergeometric functions was also given, that presents similarities with 
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the one found for the corresponding spherical modes. However, it is the 
integral form in which these mode functions are more useful, since they are 
expressed in terms of more suitable functions, and one can use them for 
orthnormalization, show completion or compute transition coefficients with 
ease. The evaluated transition coefficients between both planar and spherical 
waves, in momentum and energy bases show that there is no positive-negative 
frequency mixing between any of the discussed mode functions. 
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